Comment: Essence of intrinsic tunneling: Distinguishing intrinsic features from 

artifacts [Phys.Rev.B 72 (2005) 094503]. 
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In a recent paper V.Zavaritsky has argued that interlayer (c— axis) current-voltage characteristics 
of high temperature superconductors (HTSC) are Ohmic, and has claimed to disprove all findings 
of intrinsic tunnelling spectroscopy, as well as existence of the interlayer tunnelling and the intrinsic 
Josephson effect in HTSC, as such. 

In this comment I demonstrate, that the genuine interlayer current-voltage characteristics are 
strongly non-Ohmic, which undermines the basic postulate and the logical construction of the crit- 
icized paper. 

PACS numbers: 74.50.+r 74.72.-h 74.25. Fy 44.10.+i 



Mobile charge carriers in strongly anisotropic 
high temperature superconductors (HTSC), such as 
Bi 2 Sr 2 CaCu 2 8+( 5 (Bi-2212), are confined in Cu0 2 
planes. The out-of-plane c— axis transport in those 
compounds is caused by interlayer tunnelling which 
results in non-metallic behavior H H Q and leads to 
appearance at T < Tc of the "intrinsic" Josephson 
effect between neighboring Cu0 2 planes. At present 
all major fingerprints of the intrinsic Josephson effect 
were observed, including Fiske H II] and Shapiro 
H @ steps in Current- Voltage characteristics (IVC's); 
the Josephson plasma resonance [l(i|: thermal acti- 
vation |l ll| and macroscopic quantum tunnelling [l2T | 
from the Josephson washboard potential; and the flux 
quantization HUH 13 E3- The latter experiments 
explicitly confirmed the correspondence between the 
stacking periodicity of intrinsic Josephson junctions 
(IJJ's) and the crystallographic unit cell of Bi-2212. 
The interlayer tunnelling was also successfully employed 
for intrinsic tunnelling spectroscopy, which is unique in 
it's ability to probe bulk phonon [15| and quasiparticle 
[3 13 El 13 13 13 E3 spectra of HTSC. 

However, in a recent paper [3] and a series of other 
publications (3| V.Zavaritsky denies the existence of in- 
terlayer tunnelling in HTSC and speculates that all the 
non-linear features in c— axis IVC's are artifacts of self- 
heating. 

The two basic postulates of Refs.j3l3] a re: 

1. That the intrinsic Josephson effect does not exist, 

2. That the genuine, self-heating free c— axis IVC's of 
HTSC are Ohmic. 

The latter statement together with the assumption of 
semiconducting T— dependence of the Ohmic resistance 
leads in the presence of self- heating to non- linear IVC's, 
which according to Ref. [3] provides both qualitative 
and quantitative description of all key findings of inter- 
layer transport experiments in HTSC, proves that all the 
non-linear features of interlayer IVC's are artifacts of self- 
heating and disproves all findings of intrinsic tunnelling 



spectroscopy. 

The incorrectness of the first postulate in the face of 
overwhelming and unambiguous experimental evidence 
IlllH&UdGIlElEllUm does not require extra 
comments, except that it demonstrates an obvious failure 
of the reviewing process in several scientific journals. 

In this comment I demonstrate incorrectness of the sec- 
ond postulate. It is shown that the genuine interlayer 
IVC's of Bi-2212 are strongly non-Ohmic and are repre- 
sented by perfectly periodic, multibranch IVC's, which 
are not affected by self-heating. The huge, two order of 
magnitude, genuine non-linearity of IVC's represents the 
extent of irrelevance of the model advocated by Zavar- 
itsky to the essence of intrinsic tunnelling in HTSC and 
undermines the logical construction of Ref. [3] • 

Intrinsic tunnelling characteristics of HTSC mesa 
structures exhibit a characteristic multi-branch structure 
due to one-by one switching of IJJ's from the supercon- 
ducting into the resistive state Each time a new junc- 
tion is switched into the resistive state the dissipation 
power within the mesa increases and the effective tem- 
perature of the mesa rises as a result of self- heating [24j : 



AT = PR th , 



(1) 



where P = VI is the total dissipation power and R t h 
is the effective thermal resistance of the mesa, which de- 
pend on the sample geometry and temperature (3 13| • 
The progressive self-heating with the branch number may 
result in a systematic distortion of the IVC in the way, 
shown in Fig. 1. Namely, the critical current and the 
separation between branches decrease with the branch 
number. At large P, significant self-heating is indicated 
by progressive back-bending of the branches. 

Fig. 1 shows how the intrinsic IVC is distorted by 
extreme self- heating. It shows an IVC at the base tem- 
perature To = 4.2A" for a large Pb-doped Bi-2212 mesa 
structure, containing N ~ 40 IJJ's. The IVC is mea- 
sured in the four-probe configuration. The extremely 
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FIG. 1: ( color ) Superimposed positive and negative parts 
of the IVC of a large Pb-doped Bi-2212 mesa at T = A.2K. 
It demonstrates how the IVC is distorted by self-heating: the 
critical current and the separation between branches decrease 
with the branch number due to progressive self-heating. At 
very large dissipation power significant self-heating is indi- 
cated by progressive back-bending of the branches. 
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large P > 5mW at the end of the multi-branch structure 
in the IVC is caused first of all by a very large critical 
current density J c ~ 10 4 A/cm 2 in this Pb-doped Bi-2212 
crystal. A Similar distortion of intrinsic IVC's at com- 
parable P can be seen in Fig. 10 a) of Ref.[l5j|. Fig. 1 of 
Ref.0 

Fig. 2 shows an opposite example of a small self- 
heating. In Fig. 2 a) the IVC of a small underdoped 
Bi-2212 mesa with a small J c < 400 A /cm 2 [l<| and ap- 
proximately the same number of IJJ's, N = 34, is shown. 
Properties of this mesa were studied in Refs. |llll25l |. A 
combination of small area and J c results in a two order 
of magnitude smaller P at the end of the multi-branch 
structure, point A in Fig. 2 a), than at the similar point 
in Fig. 1. Thin lines in Fig. 2 a) represent the multi- 
ple integer of the last branch divided by N = 34, which 
indicates perfect periodicity of the branches , ty pical for 
interlayer IVC's of Bi-2212 mesas HE1HEB- Fig ' 2 
b) shows voltage jumps, Ay, between the branches as a 
function of the branch number. It is seen that AV is in- 
dependent of the branch number. Small deviations from 
the mean AV are caused by thermal fluctuations of the 
switching current fill fl^ | - 

The separation between branches provides an indepen- 
dent test for the extent of self-heating, because AV de- 
pends on T, as shown in Fig. 2 c). From Fig. 2 c) 
it follows that branches can remain periodic only if the 
successive increase of T with branch number does not 
exceed ~ 15K. This is consistent with in-situ measured 
T = 10.8K at P = 55.6/iW point A in the IVC from Fig. 
2 a) [2^|. Similarly, for the case of extreme self- heating 
shown in Fig. 1, we may conclude that the mesa is heated 
to T ~ 60K at the end of the multi-branch structure 



FIG. 2: (color) a) The IVC of a small underdoped Bi-2212 
mesa at To = 5.6K. Thin lines are multiple integers of the 
last branch divided by TV = 34 and demonstrate perfect pe- 
riodicity of the branches. Panels b) and c) show the voltage 
jumps between branches as a function of the branch number 
and the base temperature, respectively. 

at P ~ 5mW, again consistent with the measured self- 
heating AT/P ~ lOK/mW at T = 60K 25]. Further- 
more, atomic separation between IJJ's in the mesa leads 
to uniform self-heating within the mesa, see Appendix-A. 

Appendix-B presents a detailed comparison of intrinsic 
IVC's with the self-heating model of Ref. [22|. which, in 
contrast to claims of Ref. |22|, demonstrates a severe 
contradiction between experimental data and the model. 
On the other hand, this failure is obvious already from 
the raw data in Figs. 1,2. This data speaks for itself: 
The distorted periodicity of branches in Fig. 1 indicates 
significant self-heating. To the contrary, the undistorted 
perfect periodicity of branches in Fig. 2 a) implies that 
the shape of the branches is not affected by self-heating. 
However, V/I changes by approximately two orders of 
magnitude in this voltage range. 

Thus, observation of the periodic strongly non-Ohmic 
branches in interlayer IVC's can only mean that: 

i) the self-heating along the branches is negligible; 

ii) those periodic branches represent the genuine self- 
heating free IVC's of Bi-2212; 

iii) The genuine interlayer IVC's are strongly non- 
linear. 

The latter statement is the main conclusion of this 
comment. It disproves the basic postulate of Ref.j^ 
and undermines the logical construction of that paper. 
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From Fig. 2 and a large collection of similar data re- 
ported in literature, it can be concluded that the gen- 
uine non-linearity of intrinsic IVC's exceeds two orders 
of magnitude at low T, which also represents the degree 
of irrelevance of the model suggested in the criticized pa- 
per to the essence of intrinsic tunnelling in HTSC. 



APPENDIX A 

The author of Ref. speculates about "peculiar 

temperature distribution along the sample". However, 
T— distribution within the mesa can be easily estimated, 
without making any assumptions about heat flow mech- 
anism outside the mesa, or dissipation mechanism inside 
the mesa. 

Fig. 3 shows the schematics of heat flow in Bi-2212 
mesa structure. The heat P = VI, produced in the mesa, 
can flow down to the crystal and upwards into the contact 
electrode, characterized by heat resistances R cr and R e i , 
respectively. represent heat resistances of top and 

bottom parts of the mesa. The total self-heating is 



AT = PR th = P 



(Rl + i? e ;)(i?2 + Rcr) 



Rl + Rl + Re 



Rc 



(2) 



where Rth is the effective thermal resistance of the sam- 
ple, which according to Ref. is Rth = {Ah)^ 1 ~ 
4 x 10 4 K/W. The one-dimensional nature of heat flow 
within the mesa allows a straightforward estimation: 



Ri,2 — 



sN- 



1.2 



Ak c 



(3) 



where s = 15. 5A is the interlayer spacing, N±,2 is the 
number of layers in top/bottom parts of the mesa, and 
k c is the c— axis thermal conductivity. For Bi-2212, k c ~ 
Q.SW/Km at T = 30K For the mesa with A 

I i i x 10/im 2 , the thermal resistance per layer Ri(N = 
1) ~ 31K/W is three orders of magnitude smaller than 
Rth- In two extreme cases, R cr 3> R e i and R cr = R e i, 
the ratios of temperature difference per junction in the 
mesa to the total self- heating of the mesa are Ri/R t h — 
7.8 x 1CT 4 and Ri/2R th ~ 3.9 x 1CT 4 , respectively. 

For the mesa with N = 10 and P = lmW (note that 
P < 0.06mW at point A in the IVC of Fig. 2 a)) the max- 
imum temperature difference within the mesa, AT mesa is 
only - Q.3K. 

Furthermore, this is a strongly overestimated value be- 
cause: (i) k c measured in Ref.|2j| was limited by stacking 
faults in large Bi-2212 crystals, which are absent in our 
mesas. The pure phononic thermal conductivity is ex- 
pected to be almost isotropic, which would imply that 
the actual k c in the mesa is close to n a b, i.e., eight times 
higher, (ii) The estimation was obtained for the case of 
heat diffusion within the mesa. In reality, the c-axis heat 
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FIG. 3: (color) Schematics of heat flow from the mesa struc- 
ture. 



transport in Bi-2212 is predominantly phononic and bal- 
listic at the atomic scale (25[ [2t| ■ This will considerably 
reduce AT mesa for mesas with the total hight less then 
the phononic mean free path, i.e. for N < 1000. 

Therefore, for typical mesas all layers within the mesa 
are heated uniformly, irrespective of where exactly within 
the mesa the dissipation takes place. 



APPENDIX B 

One of the critical assumptions of Ref.[j^| is that 
the semiconducting zero-bias resistance Rq(T) represents 
V/I at finite bias along the IVC's. This a direct con- 
sequence of the second postulate, see above. A seeming 
matching between decaying V/I{P) and Rq(T) is demon- 
strated in Fig.l of Ref. 22] for some IVC, which, how- 
ever, does not bear any resemblance with classical multi- 
branch interlayer characteristics 0, fill Hsl flfij ] . Below I 
will repeat the same fitting for the case of conventional 
intrinsic IVC's. 

Fig. 4 represents V/I(P) (left and top axes) for the 
IVC from Fig. 2 a). The IVC in full scale is shown in inset 
of Fig. 4. The decaying part of the curve represents the 
resistance of the last branch, while the linearly increasing 
part represents multiple branches in the IVC. The dashed 
and solid lines in Fig. 4 represent Rq(T) measured with a 
small ac-modulation current and the quasiparticle resis- 
tance with all N — 34 IJJ's in the resistive state Rq (T) 
obtained from numerical differentiation of the IVC's, re- 
spectively (left and bottom axes). On the first glance, 
a good match between V/I(P) and Rq(T), in a limited 
T— range, can be obtained if we allow an arbitrary off- 
set and stretching along the T— axis of Rq(T). However, 
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FIG. 4: (color) The main panel shows the IVC from the inset, 
for the same mesa as in Fig. 2, replotted as V/I vs. P = VI 
(top axis). Dashed and solid lines represent the zero bias re- 
sistance Ro and the quasiparticle resistance Rq 4 vs. T (bot- 
tom axis), respectively. Ro(T) is matched to V/I(P) using 
the procedure suggested in Ref. 22] . Inset shows the original 
I — V characteristics of a Bi-2212 mesa at To = 5.6K. Major 
deviations between V/I(P) and Ro(T) at low P and a huge 
offset T(P = 0) ~ 70K > T = 5.6ff, indicated by an ar- 
row, reveal severe contradiction between experimental data 
and the self-heating model of Ref. [22^1 . 



several inconsistencies can be seen in such a "fit": 

i) T— offset, required for matching, leads to T(P = 
0) ~ 70K To = 5.6K, as indicated by the arrow in 
Fig. 4. A similar offset can be seen in Fig. 1 of Ref. [22| . 
Therefore, such a "fit" simply does not make any sense. 
The seeming coincidence in Fig. 4 indicates simply that 
two smooth decaying functions can be matched in the 
limited interval using two fitting parameters. 

ii) major deviations of V/I(P) and Po(T) occur at low 
P, where the main, almost 2-order of magnitude, drop of 
V/I occurs. The severe discrepancy between V/I{P) and 
Rq(T) at low bias indicates that the IVC is non- linear in 
this bias range. 

Strictly speaking, the scaling between V/I(P) and 
i?o(T) should not be expected even within the self- 
heating model because the relation between T and P in 
Eq.Q is non-linear due to T— dependence of Rth HH- 
One should instead expect self-scaling of V/I(P) curves 
at different To, using the offset AP as the only fitting 
parameter, representing the power required to compen- 
sate the difference in To between different measurements. 
Such scaling does not require constant Rth and just im- 
plies that equal P would heat the mesa by the same AT, 
irrespective of how the mesa reached the initial temper- 
ature. 

Fig. 5 shows the check for such a self-scaling for the 
optimally doped mesa (T c = 93K, A = 3.5 x 7.5/im 2 , 



a 



10 




1 P(mW) 2 



FIG. 5: (color) Check for a self-scaling of V/I(P) curves: a) a 
set of V/I(P — VI) characteristics for another Bi-2212 mesa 
at different To. The offset AP was used to fit the origin of 
each curve at P — to the nearest curve with lower To. b) 
Enlarged high P part of a), c) Temperature dependencies 
of Rq (circles) and high bias resistance R„ (squares) for the 
same mesa. It is seen that R n is smaller than Ro at any T. 



N = 10), studied in Ref. 16]. Fig. 5 a) shows a set 
of V/I(P) curves at six T from 4.9 to 220 K. For each 
curve the offset AP was chosen in such a way that the 
initial point V/I(P = 0) would coincide with V/I(AP) 
of the nearest curve with lower Tq. Apparently, the self- 
scaling is non existing for conventional intrinsic IVC's. 
Fig. 5 b) shows zoom-in of the high power part of the 
plot. Here self-heating should be the largest. It is seen 
that the IVC's asy mptotically approach some "normal" 
resistance R n [16J. Experimental values of R n and Po 
vs. To are shown in Fig. 5 c). Comparison of R n and 
Po shows that P„ is smaller than Po at any T. This 
makes the fitting of V/I by Po(T) technically impossible, 
indicating that the IVC's have to be non-Ohmic. 
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